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We investigate the size of potential higher pseudoscalar resonance contaminations of the esti- 
mates of isospin-conserving and isospin-violating ttNN couplings obtained using the 3-point function 
method in QCD sum rules. For the isospin-conserving case it is shown that conventional models 
of the isovector pseudoscalar spectral function imply resonance decay constants large enough to 
create significant contaminations, and that assuming these models are incorrect, and that the decay 
constants are actually much smaller, implies physically implausible values for the flavor-breaking 
quark condensate ratios. For the isospin-violating case it is shown explicitly that such resonance 
contamination is unavoidably present and precludes using the 3-point function method as a means 
of estimating the at present unmeasured isospin-violating ttNN couplings. 
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I. INTRODUCTION 



The general framework of QCD sum rules [1-3] has recently proven popular as an approach to the computation of 
observable of relevance to problems in nuclear and few-body physics which has, in contrast to many effective hadronic 
models, a rather direct connection with QCD itself. The method is attractive, first, because it is based on rather 
general properties of the underlying field theory (the operator product expansion (OPE), analyticity, unitarity, and 
the existence of appropriately subtracted dispersion relations) and, second, because it provides a means of relating 
integrals over physical spectral densities to the behavior of correlators at large spacelike momenta (obtained via the 
OPE) which, via the Borel transformation of the original dispersion relation, simultaneously exponentially suppresses 
higher energy contributions to the physical spectral integrals and factorially suppresses the contributions to the 
OPE associated with higher dimension operators. This means that one has the quite reasonable hope of, in favorable 
circumstances, constructing a sum rule relating observable parameters (such as masses, couplings and decay constants) 
occuring on the phenomenological side of the sum rule to a small number of vacuum condensates (which parametrize 
non-perturbative effects in QCD) appearing on the OPE side of the sum rule. 

In the case of 2-point functions, it is rather easy to see what "circumstances" are favorable to such an analysis. First, 
one should know that one is considering a channel where the relevant spectral function consists of a lowest lying single 
resonance contribution well-separated from higher resonance and/or continuum pieces. This allows one to choose a 
Borel mass to strongly suppress the contributions to the weighted spectral integral from the more complicated part 
of the physical spectral function, and hence obtain an expression for the phenomenological side of the sum rule that 
is dominated by a few observable parameters associated with the lowest lying physical state. Second, for such Borel 
masses, it must simultaneously be the case that the OPE side of the sum rule is well-converged at operators of low 
enough dimension that the corresponding vacuum values are already known from other analyses. Leinweber [4] has 
provided a very clear discussion of the criteria to be satisfied for the applicability of the sum rule method, and also a 
procedure for sensibly checking the validity of these criteria in a given case. 

In the case of observables like the isospin-conserving and isospin-breaking nNN couplings, the situation is slightly 
more complicated. In the past, two different approaches have been taken to estimating these couplings in the QCD 
sum rule framework. 

The first of these approaches is the 3-point function method [2,5-7]. In this approach one considers the 3-point 
vacuum correlator 

G-^^NNiPi:P2,q) = J d^xid'^X2exp {ipixi ~ ip2X2) {0\T {xNixi)Pj=iiO)x'Nix2))\0) , (1) 

where a labels the isospin of the pion (± or 0), iV stands for either proton or neutron, Pj^i = is a pseudoscalar 
interpolating field for the pion, and xn is the loffe interpolating field for the nucleon. In what follows we will use the 
following notation for the pseudoscalar currents: 
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Pf = qfil^qf 

P+ = V2uij5d, (2) 

where / is a flavor label f = u,d, s. The a = component of the isovector current multiplet is, as usual, 

Pi = ui'y^u — di'y^d . (3) 

In what follows we will also require the flavor neutral isoscalar currents. For these we can use either the strange and 
light quark combinations 

Pq = ui^^u + di'y^d 

Pq = si-y^s (4) 

or the singlet and octet combinations 

Pq = [uij^u + dij^d + sij^s] / Vs 

Pq = [uij5U + dijsd - 2si^5s\ /Vs (5) 

as a basis. The loffe currents for the nucleons are deflned by 

Xp = eabc [{KfC^y) jr,rd'] , (6) 
Xn = eabc [{[d-fCj^d") ^.ru"] . (7) 

Schematically, one then analyses the correlator using both a hadronic model (which involves the relevant t:NN 
couplings as parameters) and the OPE, and matches the two representations, appropriately transformed, to extract 
the wNN couplings. In order to perform this analysis, it is necessary that the momentum q on the external pion 
leg of the correlator be large and spacclikc, otherwise the OPE of the correlator will not be valid truncated to the 
low-dimension operators for which the analysis is practical. This condition, however, means that one is rather far from 
the pion pole. This problem is dealt with by looking for terms in the OPE which have the same Lorentz structure as 
the pion contribution and in addition have a pole of the form 1/q'^. However, as stressed by Birse and Krippa [8], this 
procccdurc is inherently rather dangerous, since it is not a priori possible to separate contributions from the pion and 
from higher resonances to the coefficient of l/q^. Since, moreover, one must certainly work at = —q^ > 1 GeV^, 
i.e., rather far from the pion pole, one may no longer reasonably count on the proximity of this pole to conclude 
that the pion contribution is dominant, as it would be for small Q^. The reliability of the 3-point function method 
thus rests crucially on the assumption that the higher resonance contributions, in the region of values under 
consideration, are small. The plausibility of this assumption has not previously, to our knowledge, been investigated. 

The second approach to the ttA^A'^ couplings is the 2-point function method [2,6,8,9]. In this approach, one considers 
the 2-point correlator 

n(p,g) = ij d^x exp(ipx)(0|T(xjv(x)xiv(0)) {n^q)) . (8) 

For large spacelike values of p, the OPE for the product of the two nucleon interpolating flelds can presumably be 
truncated at operators of relatively low dimension. One may then look at the vacuum-to-pion matrix elements of 
these operators in order to estimate the ttNN couplings. As has been stressed by many authors, in order to be 
able to remove contributions associated with N N* transitions, one should look not at the 75 term in the OPE, 
but rather the structure ^75 [8,11-14]. The advantage of the 2-point method is, of course, that it completely avoids 
the problem of potential contamination from higher resonance contributions, which is unavoidable in the 3-point 
function method. The disadvantage of the method is that, although the vacuum-to-pion matrix elements of the lowest 
dimension operators can all be evaluated with good accuracy using chiral perturbation theory (ChPT) [10], those for 
the higher dimension operators are less certain. Thus, for example, the uncertainty in the value of the matrix element 

g,(0|JG^'^7.'ik+(<z)) (9) 

leads to uncertainties of ~ 15% in the prediction for the isospin-conserving ttNN coupling using the 2-point function 
method [8]. 

To date, the isospin-conserving ttNN coupling has been considered using both the 2-point [2,6,8,9] and 3-point 
[2,5,6] function methods, and the isospin- violating ttNN couplings using the 3-point function method [7]. The isospin- 
conserving coupling is, of course, known experimentally, so the point of computing it using QCD sum rules is primarily 
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to test the plausibility of the assumptions and truncations that go into the evaluation. The hope is that success in 
computing the coupling in the isospin-conserving case might serve as an indication of the reliability of the approach 
employed and hence make the analogous calculation of the unmeasured isospin-violating couplings also plausibly 
reliable. Since both the 2-point and 3-point function treatments are successful in this regard, albeit it with significant 
theoretical errors, it would seem reasonable to attempt to proceed to the isospin-violating case using either method. 
A first attempt in this direction was made by Meissner and Henley [7] , employing the 3-point function method. Since 
the isospin violating coupling has not yet been experimentally determined, this estimate of is considerable potential 
interest, particularly in view of the recent revival of interest in the question of isospin-breaking in few-body systems 
(see for example the review of Rcf. [15] for an extensive discussion of the situation up to 1990, and Ref. [16] for a list 
of more recent papers on the subject). We would, therefore, like to understand whether, given the potential problems 
of the 3-point function method, this estimate is reliable or not. 

In the present paper, therefore, we will investigate the question of higher resonance contaminations in the 3-point 
function method, which is the biggest potential roadblock to using the method in the isospin-violating case. We 
will show that existing (albeit model-dependent) understanding of the spectral density in the isovector pseudoscalar 
channel implies that significant contamination is present, even in the isospin-conserving case, and that requiring this 
understanding to be incorrect, and the contamination to be small, is equivalent to rather strong (and physically 
implausible) statements about the values of flavor-breaking ratios of quark condensates. We will then proceed to 
show that certain features of the isospin-violating analysis itself also clearly indicate the presence of significant higher 
resonance contamination, implying that the 3-point function method cannot be reliably employed to extract the 
isospin-violating ttA^'A'' couplings. 



II. THE ISOSPIN-CONSERVING ANALYSIS 



One does not know, a priori, the size of the couplings of the excited isovector pseudoscalar mesons, 7r(1300), 
7r(1800), • • • to the nucleon. Thus, in order to be certain that the 3-point hmction method for extracting g„NN is not 
contaminated by contributions from these resonances it is necessary that the conditions 

«77#^ (10) 



be satisfied, where /m is the decay constant for meson M and tt' stands for any of the excited 1=1 pseudoscalar 
mesons. 

The excited state pseudoscalar decay constants are not known experimentally and, owing to the fact that they have 
a chiral suppression, [1] will be very difficult to measure. (The 7r(1300) decay constant could, in principle, be measured 
by separating the small pseudoscalar contribution from the large overlapping ai contribution in r decays via a detailed 
spin-parity analysis [17].) However, they are related to the light current quark mass combination mu + md via a series 
of finite energy sum rules [18]. The most recent analysis of this mass combination [18] models the continuum part 
of the spectral function in the isovector pseudoscalar channel in terms of a sum of 7r(1300) and 7r(1800) resonances. 
The relative contribution of the two resonances in the model continuum spectral function is constrained by duality, 
and the overall normalization is set by assuming the threshold value of the continuum spectral function, which can be 
estimated using ChPT at tree level, is saturated by the tails of the resonance contributions. If we assume the model 
spectral functions so constructed (which are tied to the usually quoted values of the MS masses) are reasonable, then 
we can read off the corresponding values of the decay constants for the 7r(1300) and 7r(1800). Taking the model with 
the best duality fit from Ref. [18], we find 

/^(i30o) = 2.2 MeV 

/.(1800) = 1-0 MeV. (11) 

Forming the products of the couplings, fMinl^, to the pseudoscalar 7=1 current, and the propagators, l/{Q'^ + m%[), 
evaluated at ~ 1 GeV^, we then find 



= 1.8, 2.7, and 0.76 MeV (12) 

Q2=i GeV^ 



for M = TT, 7r(1300), and 7r(1800), respectively. For such values of the excited pseudoscalar decay constants, therefore, 
one would be forced to conclude that the contamination from, certainly the 7r(1300), and most likely also the 7r(1800), 
would be far too large to make the method reliable. 
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Of course, one might object that the above argument, relying as it does on the model spectral functions of Rcf. [18], 
is model-dependent and therefore not conclusive. Indeed, the validity of the method of Ref. [18] for setting the overall 
normalization (and hence the overall scale of the corresponding decay constants in Eq. (11)) has been questioned [19], 
leading to the suggestion that the normalization might actually be significantly smaller than employed in Ref. [18]. 
However, even if the normalization of the continuum spectral function were to be decreased by an order of magnitude, 
the corresponding decay constants would be decreased only by a factor ~ 3, leaving the product of the 7r(1300) 
coupling to the current and propagator at the level of ~ 50% of that of the tt. 

It thus appears extemely unlikely that the couplings of the higher pseudoscalar resonances can be neglected in the 
3-point hmction analysis of g^rWAr- Since, however, one does not actually have an experimental value for /7r(i3oo)! one is 
still free to imagine that the decay constants are, for some reason, much more than just an order of magnitude smaller 
than those corresponding to the model spectral functions of Ref. [18] (for example, a spectral function a factor of 60 
times smaller than that of Rcf. [18] would bring the 7r(1300) coupling times propagator factor to below 20% of the 
corresponding product for the tt). This appears a rather unlikely prospect, but is one that cannot be presently ruled 
out. However, it is important to realize that such an assumption is not is not without other non-trivial consequences. 

Let us, therefore, for the moment, accept the (albeit unlikely) prospect of extremely small excited pseudoscalar 
decay constants and consider the consequences of such an assumption for the quark masses and condensates. The 
first consequence is obvious, namely, if we strongly suppress the continuum contribution to the finite energy sum rule 
analysis for the quark mass, then we commit ourselves to significantly lower values of rriu + rud, of order 6 MeV (in 
the MS scheme, at a scale of 1 GeV^, this value to be compared with the conventionally quoted value 12 MeV). This 
is not necessarily a problem, since a recent analysis of world lattice data also suggests a significantly lowered value of 
m„ + nid [20], a possibility also noted in Ref. [19]. It does, however, force one to significantly larger values of the light 
quark condensate, which can cause problems for the stability of the sum rule for the nucleon mass [21]. Moreover, 
such an assumption actually corresponds to rather strong constraints on the ChPT low energy constant (LEG), H2, 
(where we adhere to the notation of Gasser and Leutwyler [10] throughout), which LEG governs the fiavor breaking of 
the quark condensates [10]. The reason for the existence of this constraint is that the inverse weighted (corresponding 
to n = — 1 in the notation of Ref. [18]) finite energy sum rule for the correlator 

*5(g2) = if d'^xe'i-''{0\T{d^'Al-\x),d''Ai+\0)}\0), 

= {md + mufij d4a;e'«-^(0]T{P(-)(a;),P(+)(0)}]0) (13) 

with A^f^^ the charged isovector axial current and p(+) the corresponding charged isovector pseudoscalar current, can 
be re- written as a sum rule for the continuum portion of the pseudoscalar spectral function, as follows: 

r dtlm^5{t) 3 / ^^2 n , p / M 8to^/^ 

/ = -^[mu[s)+md{s)) s[l + Ro{s)] p^i^L^-H^) (14) 

where m„^d(s) are the MS running masses at scale s, Ro{s) contains perturbative corrections [18], F is a leading order 
GhPT LEG, equal to the pion decay constant in the chiral limit, and 2Lg — iJJ is a scale-independent combination 
of fourth order GhPT LEG's. If we suppress the integral on the LHS by a large factor like 60, and also the running 
masses by of order a factor of 2, it turns out that we drive 2Lg — H2 to values more than 3 times smaller than those 
obtained in Ref. [18]. This in turn implies that i?2 ("^r?) must necessarily be positive. We now argue that such values 
for H2{m^) lead to physically implausible predictions for the ratios of quark condensates. 

To see this, note that, once H2 is fixed, the flavor breaking ratios of quark condensates are simultaneously fixed at 
next-to-leading order in the chiral expansion [10], for example, 

^Q|^^|Q^ - 1 -h 3/i„ - lUK - ^lr, + (^^8 + ^2) , (15) 

where /iM = ml^j \og{m\f / fi'^) / 32Tr^ , with ^ the GhPT renormalization scale. With H2{mn) > and Lg(m,,) = 
(1.1 ± 0.3) X 10-3 [22], we see from Eq. (15) that 

j» > 1.38 . (16) 

Thus, assuming that the excited resonance decay constants are sufficiently small to be able to neglect their contri- 
butions to the 3-point function sum rule simultaneously commits one to the highly unnatural situation of a strange 
quark condensate larger in magnitude than the light quark condensate. In addition, one finds that, owing to the 



4 



relation between the flavor breaking and isospin breaking condensate ratios given by Eq. (9.5) of Ref. [10], such a 
value for the strange to up quark condensate ratio implies, for the isospin-breaking condensate ratio, 7, defined by 

_ {0\dd\0) - {0\uu\0) 
^~ {0\uu\0) ' ^ ' 

a value 

7 > 1.5 X 10-3 , (18) 

in contradiction with extractions of 7 from a variety of sources, [3,23-27] all of which obtain 7 < 0. 

In view of the results of the last paragraph, we conclude that the hypothesis that one may neglect the higher reso- 
nance contamination in the 3-point function analysis of the isospin conserving coupling Qt^nn is a highly unpalatable 
one. In the next section, we will turn to the case of the analysis of the isospin-breaking couplings and demonstrate 
more directly the presence of analogous higher resonance contaminations. 



III. THE ISOSPIN- VIOLATING ANALYSIS 



In this section we will concentrate on the 3-point function analysis of the difference of n^pp and n^nn couplings, 

^9 = 9-KOnn - 9%0pp- (19) 

In order to perform this analysis it is necessary to take into account the fact that Pi is not a suitable interpolating 
field for the physical 7r°, if one wishes to treat isospin breaking effects. This follows from the observation that 



(20) 



As a consequence, if one were to use Pj* as tt'^ interpolating field, then even if one could ignore higher resonance 
contributions, the result of the analysis would be a mixture of the isospin breaking tt coupling and the isospin 
conserving r] coupling (the latter multiplied by an isospin breaking factor describing the coupling of the t] to the 
J = 1, = current). As noted above, there is no means to separate the contributions corresponding to different 
mesons in the 3-point function approach. Meissner and Henley [7] , who first performed the isospin violating analysis, 
dealt with this problem by choosing a current combination with no vacuum-to-7? matrix element, namely, 



where the pseudoscalar currents are as defined above and the choice 

€ = 60 = — — 

4 rris — m 

where to = (to„ -|- md)/2, and is the leading order tt-t] mixing angle [10], ensures that 

(0|P,o|,7) =0 



(21) 



(22) 



(23) 



to leading order in the chiral expansion. The choice of interpolating field with this property is not unique [16]; in fact, 
for any a, defining 



P(a) = (P„ - Pd) + e(a) a(P, 

1 



(P„ - Pd) + e(a) 



V3 



Pd) + (l+a)P, 
-Ps + {3a + l)Po 



one may find an e{a) such that 



< 0]P(a)|7? >= . 



(24) 



(25) 



The set of solutions of Eq. (25), as a function of a, have been worked out to next-to-leading order in the chiral 
expansion in Ref. [16]. Among other results of this analysis, it is found that the next- to- leading corrections are 
significant; for example, at next-to-leading order, the Meissner-Henley field choice must be modified to 
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P^o + 1.27^0^0 (26) 

if one wishes to maintain zero vacuum-to-r; matrix clement. 

We will now explain why the existence of the above family of potential tt'^ interpolating fields is of relevance to 
our current discussion. Note that the choice of Pt^o was made by Meissner and Henley with the problem of potential 
higher resonance contamination in mind. Indeed, they argued that this choice of interpolating field is the one that 
would suppress possible rj' contributions [7]. It turns out, however, that this is not the case, and the reason that it 
fails to be so leads us immediately into a consideration of the larger class of 7r° interpolating fields. 

Let us, therefore, first outline why the interpolating field P^o necessarily has non- vanishing coupling to the rj' . In 
the chiral limit, of course, there is no flavor breaking whatsoever, hence no isospin breaking couplings, and no flavor 
breaking meson decay constants. Once we introduce the quark mass matrix, with its flavor breaking rUg — rh difference 
and isospin breaking — niu difference, all flavor and isospin breaking effects are potentially present. In the case of 
the isospin breaking and flavor breaking decay constants of the r?', one may obtain a leading estimate of the ratio of 
these decay constants using SU{i)F arguments. Indeed, we know that the breaking is produced by the quark mass 
matrix, which has the following decomposition into singlet, octet isovector and octet hypercharge components 

M =^{ms + 2m) - ^(m^ - m„)A3 - -^(mg - rh)\s, . (27) 

To leading order in the isospin-breaking and flavor-breaking mass differences, therefore, the vacuum-to-77' matrix 
elements of P° and Pg are given simply by the product of the mass-dependent coefficients of A3 and Ag in Eq. (27) 
with a common 8^x8^—^ If reduced matrix element. Recasting the ratio of mass factors in terms of the mixing 
angle ^0 deflned above, one then finds, straightforwardly, that 

<0|P„o|r/'>=36lo<0|Po^|77'> (28) 

where < 0|Pq jry' >= 0{ms — fa). The RHS is thus non-zero, of 0{rad — ruu), and in fact has a munerical enhancement 
(the factor 3) brought about by the fact that the couplings of the rj' to the P° and Pq components of P^o add 
coherently. Note that a similar argument, using a first order treatment of flavor and isospin breaking, would predict 
that the isospin violating and flavor violating axial ciirrcnt 77' decay constants were of the same sign, in agreement 
with the results of a recent QCD sum rule analysis of the isospin violating < 0|T(A^A^)|0 > correlator [28]. 

If one considers allowing an admixture of the singlet pseudoscalar current into the interpolating fleld (i.e, allowing 
a to deviate from —1/3), one has, of course, in addition to the 8^x8^— reduced matrix element which governs 
< OlP}*!?]' > and < 0|Po|?7' >, the 1^ x I/t- ^ Ip reduced matrix element relevant to < 0|Pq|?7' >. One can then 
certainly, in principle, find a value of a such that < 0\P{a)\ri' >= 0. The problem is that, even to do so at leading 
order in the quark masses, one would need to know the ratio of the two reduced matrix elements above, and this 
information is not available. Moreover, even if it were, this would not necessarily ensure that, for such a value of a, 
the couplings of the higher resonances other than the rj' were small for the same value of a. 

Since we do not know, a priori, how to choose a 7r° interpolating field (i.e., a value of a) to remove even 77' 
contamination, let alone possible contamination associated with yet higher pseudoscalar resonances, it is necessary 
to look for some sort of post facto indication of the absence of such contributions. One obvious way of doing so 
is to study the extracted results for what is nominally the isospin violating ttNN coupling, 6g, as a function of a. 
If one can find a region of a values for which the results are not sensitive to a, then one might argue this was a 
signal that, for such values of a, the effect of couplings to all higher resonances is negligible. In contrast, if one is 
unable to find such a region of a values then, since the various interpolating fields differ only in their couplings to 
the excited pseudoscalar resonances, beginning with the rj', it is clear that, in general, there are large contaminations 
from the higher pseudoscalar resonances and that, as a consequence, one has no reliable way of choosing a particular 
interpolating field for which such contaminations are small. We will see, unfortunately, that it is the latter situation 
which holds for the 3-point function analysis. Moreover, we will demonstate that the variation of 5g with a is 
essentially as large as the value of Sg extracted in the Meissner-Henley analysis, and hence that no reliable estimate 
of Sg can be made using the 3-point function method. 

In order to demonstrate the claims of the last paragraph, it is necessary to understand how the generalization from 
the specific Meissner-Henley interpolating field choice, a = —1/3, to arbitrary a affects the sum rule for the isospin 
breaking ttNN coupling. It is straightforward to show that, in the general case, the final sum rule (the analogue of 
Eq. (24) of Ref. [7]) becomes 



Sg ] 2 4 /6Mm 
= 7 + —Kia) + 

{2t:)VU^M^ 

M6 + i£,2 < > ^2 



md-mu\ , 6Mn Mn 



21 



rrid + muj Mn 



(29) 



6 



where M is the Borel mass, < > the gluon condensate, g,, the strong coupling constant, SA-In the quark-mass- 
difference contribution to the nucleon mass splitting, and K{a) the coefficient of Pu + Pd appearing in the interpolating 
field P{a) (so, for example, for the Meissner- Henley choice P^o, k = 60 /V^)- The terms in the second line of Eq. (29) 
arise from evaluating the isospin breaking difference of couplings of the loffe currents to the neutron and proton 
states using the chiral odd sum rule for the nucleon two-point functions [7]. The reason for the appearance of K{a) in 
Eq. (29), i.e., the dependence on only the light quark / = content of the interpolating field, is that contributions of 
the strange quark component of the interpolating field to the OPE side of the sum rule vanish to the order considered 
in obtaining the sum rule. 

From the expression for P{a) in Eq. (24), it is evident that 

K{a) = ae{a) . (30) 

In order to complete our investigation we, therefore, require only the expression for e(a), obtained in Ref. [16]: 



e(a) = 



1- 

32 
3^2 



-(10 + 9a)nn + 6(1 + a)l^K + (4 + Sa)^,, 
3P2 



(4 + 3a) {fnl - ml) (3L; + Lg) + (^^|^ 




1 + 



^^|^(l + log«/M^)) 



(31) 



where m|- is the average of the and squared masses, the chiral log terms /Um are as defined above and the 
are the usual renormalized fourth order LEC's, in the notation of Gasser and Leutwyler [10]. In Eq. (31) the 
expression has been written so that the 1 occurring in the square brackets corresponds to the contribution obtained at 
leading order in the chiral expansion, while the remaining terms give the next-to-leading order corrections. The above 
results reproduce the Meissner-Henley field choice for a = —1/3, if one keeps only the leading order contribution in 
Eq. (31). 

If we take the latest evaluation of the quark mass ratios from ChPT [29] , then we have 



9o = (1.1 ±0.2) X 10" 



(32) 



(The upper end of the error bound would correspond to the evaluation of r = (m^ — m„)/(TOd + '7i„) obtained using 
77 Stt, which corresponds to large values for the violation of Dashen's theorem advocated by a number of authors 
[30 35].) For the Meissner-Henley choice a = —1/3, we then find, for the contribution to —Sg/ g^^NN generated by 
the isoscalar component of the 7r° interpolating field (required to remove the r] contamination from the final result), 



[Sg/g^NN]^ 



(-1/3) 



4_^ 
3^3 



[1 + 0.27] = 1.1 X 10" 



(33) 



while, for comparison, for the choice a = — 1, which removes the strange quark content from the interpolating field, 
we find 



-Sg/g 



ttNN 



k{-1) 



-^6»o[l + 0.21] = 3.1 X 10" 
o 



(34) 



To understand the implication of these results, one should bear in mind that the range of values for —Sg/girNN 
extracted by Meissner and Henley was (1.7-3.0) x 10"^. (The range reflects the full range of uncertainties in all of 
the input parameters, 7, Oq and 5Mjq.) We thus find that the corrections required to remove the t] contamination 
are large on the scale of the result obtained. Moreover, since K{a) = ae{a), and the results of Ref. [16] show 6(0;) 
to be slowly varying with a, and greater than 1 over a wide range of a values, we see also that the results for what 
is nominally the isospin breaking coupling is very sensitive to a, varying, for example, by an amount as large as the 
maximum value quoted by Meissner and Henley over the range between a = and a = —1. We thus conclude, in light 
of the discussion above, that the 3-point function evaluation of the isospin breaking coupling is plagued by unknown 
higher resonance contamination, and as such cannot provide a reliable estimate of this quantity. 
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IV. CONCLUSIONS 



We have shown that the 3-point function method for the treatment of both the isospin conserving and isospin vio- 
lating ttNN couplings is plagued by problems with higher resonance contamination. In the course of this investigation 
we have also seen how information from the finite energy sum rules for the light quark masses, chiral perturbation 
theory, and sum rules for the chiral LEC's can sometimes be profitably employed to elucidate the physical content of 
other sum rules treatments. 
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